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1 Introduction 


If neutrinos are massive and mixed in the SM, due to the smallness of the neutrino masses, 
the predicted branching ratios for the lepton flavor violation (LFV) are so tiny that they are 
completely unobservable. On the other hand, in the supersymmetric framework the situation 
is quite different. Many authors have already studied the LFV in the minimal supersymmetric 
standard model (MSSM) with right-handed neutrinos assuming the relevant neutrino mass ma¬ 
trix iiii- In the MSSM with soft breaking terms, there exist lepton flavor violating terms 
such as off-diagonal elements of slepton mass matrices and trilinear couplings 

A|j. Strong bounds on these matrix elements come from requiring branching ratios for LFV 
processes to be below observed ratios. For the present, the most stringent bound comes from 
the /U ^ e -|- 7 decay (BR(^ —> e -|- 7) < 1.2 x 10 “^^) p. However, if the LFV occurs at tree 
level in the soft breaking terms, the branching ratio of this process exceeds the experimental 
bound considerably. Therefore one assumes that the LFV does not occur at tree level in the 
soft parameters. This is realized by taking the assumption that soft parameters such as 

Afj, are universal i.e., proportional to the unit matrix. However, even though there 
is no flavor violation at tree level, it is generated by the effect of the renormalization group 
equations (RGE’s) via neutrino Yukawa couplings. Suppose that neutrino masses are produced 
by the see-saw mechanism |^, there are the right-handed neutrinos above a scale M^. Then 
neutrinos have the Yukawa coupling matrix Yj^ with off-diagonal entries in the basis of the 
diagonal charged-lepton Yukawa couplings. The off-diagonal elements of drive off-diagonal 
ones in the and Afj matrices through the RGE’s running |^. 

One can construct Y^/ by the recent data of neutrino oscillations. Assuming that oscillations 
need only accounting for the solar and the atmospheric neutrino data, we take the LMA-MSW 
solution of the solar neutrino. Then, the lepton mixing matrix, which may be called the MNS 
matrix or the MNSP matrix [^, is given in ref.[lR|. Since the data of neutrino oscillations only 
indicate the differences of the mass square Arnfj, neutrinos have three possible mass spectra: 
the hierarchical spectrum 171^3 3 > m^2 ^ , the quasi-degenerate one m^i ~ m^2 — and 

the inverse-hierarchical one rriyi ~ 771^2 ^ 771^3. 


2 LFV in the MSSM with Right-handed Neutrinos 

The superpotential of the lepton sector is described as follows: 

hhiepton = Y^LHdel + YyLHuvl + ]^ufm^ul, (2.1) 

where Hu,Hd are chiral superfields for Higgs doublets, L is the left-handed lepton doublet, 
and Vji are the right-handed charged lepton and the neutrino superfields, respectively. The Yg 
is the Yukawa coupling matrix for the charged lepton, Mj^ is Major ana mass matrix of the 
right-handed neutrinos. We take Yg and Mj^ to be diagonal. 

It is well-known that the neutrino mass matrix is given as uiy = (YyVy)^ C^uVu) > via 
the see-saw mechanism, where Vu is the vacuum expectation value (VEV) of Higgs Hu. The 
neutrino mass matrix niy is diagonalized by a single unitary matrix = U^j^gmi^UMNS) 

where Umns is the lepton mixing matrix. Following the expression in ref. 0 , we write the 



neutrino Yukawa coupling as 


Y, = R uLs , ( 2 . 2 ) 

where R is a 3 x 3 orthogonal matrix, which depends on models. 

At first, let us take the degenerate right-handed Majorana masses M^i = = M/j. 

This assumption is reasonable for the case of the quasi-degenerate neutrino masses. Otherwise 
a big conspiracy would be needed between and Mj^. This assumption is also taken for cases 
of the inverse-hierarchical and the hierarchical neutrino masses. We also discuss later the effect 
of the deviation from the degenerate right-handed Majonara neutrino masses. Then, we get 
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It is remarked that Y))Y^, is independent of R in the case of = M^2 = 

As mentioned in the previous section, there are three possible neutrino mass spectra. The 
hierarchical type {rriyi <C my2 ^ gives the neutrino mass spectrum as 
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and the inverse-hierarchical type {mui ~ mu2 S> mus) gives 
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We take the typical values Aml^^ = 3 x 10 ^eV^ and AttIq = 7 x 10 ®eV^ in our calculation 
of the LTV. 

We take the typical mixing angles of the LMA-MSW solution such as S23 = l /\/2 and S12 = 
0.6 1^, in which the lepton mixing matrix is given in terms of the standard parametrization of 
the mixing matrix |j^. The reacter experiment of CHOOZ |]^ presented a upper bound of S13. 
We use the constraint from the two flavor analysis, which is S13 < 0.2 in our calculation. If we 


take account of the recent result of the three flavor analysis \V2^, the upper bound of S13 may 
be smaller than 0 . 2 . Then, if we use the results in our results of /r —> e -|- 7 are reduced at 
most by a factor of two. In our calculation, the CP violating phase is neglected for simplicity. 

Since SUSY is spontaneously broken at the low energy, we consider the MSSM with the soft 
SUSY breaking terms: 

-Tsoft = (m|)iAl 4 ' + (mDijeUe^j + + Af-R^eYL^ + , ( 2 . 7 ) 


where m|,m| and m? are mass-squares of the left-handed charged slepton, the right-handed 
charged slepton and the sneutrino, respectively. Ag and Aj^ are A-parameters. 











Note that the lepton flavor violating processes come from diagrams including non-zero off- 
diagonal elements of the soft parameter. In this paper we assume the m-SUGRA, therefore we 
put the assumption of universality for soft SUSY breaking terms at the unification scale: 

= (m|)y = = • • • = 5ijml , A!" = Y,,aomo, A® = Ygaomo , (2.8) 


where mo and oq stand for the universal scalar mass and the universal A-parameter, respectively. 
The RGB’s for the left-handed slepton soft mass are given by 
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while the first term in the right hand side is the normal MSSM term which has no LFV, and 
the second one is a source of the LFV through the off-diagonal elements of neutrino Yukawa 
couplings. 


3 Numerical Analyses of Branching Ratios 


Let us calculate the branching ratio of e* —> ej -|- 7 (j < i). The decay rate can be calculated 
using as 

r{ei ej + 7 ) = • ( 3 - 1 ) 

Since we know the relation m^, » mg^., then we can expect \A^\ » \A^\. The decay amplitude 
is approximated as 
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where 02 is the gauge coupling constant of SU{ 2 )i and ms is a SUSY particle mass. The RGB’s 
develop the off-diagonal elements of the slepton mass matrix and A-term: 
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where Mx is the GUT scale. 

We present a qualitative discussion on (yJy,) 21 before predicting the branching ratio 
BR(^ ^ e -|- 7). This is given in terms of neutrino masses and mixings at the electroweak 
scale as follows: 

(YtY,)2i = ^ [U^2U:2(m,2 - m,i) + U^3U:3{m,3 - m,i)] , (3.4) 

U/ 


where Vu = vsinfd with v = 174 GeV is taken as an usual notation and the unitarity condition of 
the lepton mixing matrix elements is used. Taking the three cases of the neutrino mass spectra, 



one obtains the following forms, repectively, 
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where we take the maximal mixing for the atmospheric neutrinos. Since Ue 2 — for the 

bi-maximal mixing matrix, the first terms in the square brackets of the right hand sides of 


eqs. (|3.5D can be estimated by putting the experimental data. For the case of the degenerate 
neutrino masses, (yJYi/) 2 i depends on the unknown neutrino mass scale rriu- As one takes the 
smaller m^, one predicts the larger branching ratio. In our calculation, we take = 0.3eV , 


which is close to the upper bound from the neutrinoless double beta decay experiment |W], and 
also leads to the smallest branching ratio. 

We also note that the degenerate case gives the smallest branching ratio BR(/i —>6 + 7 ) 
among the three cases as seen in eqs.(^) owing to the scale of rriu- It is easy to see the fact 
that the second terms in eqs. (|3.5[) are dominant as far as ~ 0.01 (degenerate), 0.01 (inverse) 
and 0.07(hierarchy), respectively. The magnitude and the phase of are important in the 
comparison between cases of the inverse-hierarchical and the normal hierarchical masses. In the 
limit of t/e 3 = 0 , the predicted branching ratio in the case of the normal hierarchical masses is 
larger than the other one. However, for t/es — 0.2 the predicted branching ratios are almost the 
same in both cases. 

At first, we present numerical results in the case of the degenerate neutrino masses assuming 
Mr = Mfll. The magnitude of M^ is constrained considerably if we impose the b — T unfication 
of Yukawa couplings In the case of tan/3 < 30, the lower bound of Mr is approximately 


lO^^GeV. We take also Mr < lO^^GeV, in order that neutrino Yukawa couplings remain below 
0(1). Therefore, we use Mr = 10^^, lO^^GeV in our following calculation. 

We take a universal scalar mass (m-o) for all scalars and oq = 0 as a universal A-term at the 
GUT scale {Mx = 2 x 10^® GeV). The branching ratio of /r ^ e-l -7 is given versus the left-handed 
selectron mass rrie^ for each tan /3 = 3, 10, 30 and a fixed wino mass M 2 at the electroweak 
scale. In fig.|I|, the branching ratios are shown for M 2 = 150, 300 GeV in the case of Ues = 0.2 
with Mr = lO^'^GeV and = 0.3eV, in which the solid curves correspond to M 2 = 150GeV 
and the dashed ones to M 2 = 300GeV. The threshold of the selectron mass is determined by 


the recent LEP2 data IQ for M 2 = 150GeV, however, for M 2 = 300GeV, determined by the 
constraint that the left-handed slepton should be heavier than the neutralinos. As the tan f3 
increases, the branching ratio increases because the decay amplitude from the SUSY diagrams 
is approximately proportional to tan /3 |j|] . It is found that the branching ratio is almost larger 
than the experimental upper bound in the case of M 2 = 150GeV. On the other hand, the 
predicted values are smaller than the experimental bound except for tan/3 = 30 in the case of 
M 2 = 300GeV. 

Our predictions depend on Mr strongly, because the magnitude of the neutrino Yukawa 





















coupling is determined by as seen in eq.(^). If reduces to lO^^GeV, the branching 
ratio becomes 10^ times smaller since it is proportional to MH. The numerical result is shown 
in fig|. 

Next we show results in the case of the inverse-hierarchical neutrino masses. As expected in 
eq.(^]^), the branching ratio is much larger than the one in the degenerate case. In fig.|3|, the 
branching ratio is shown for M 2 = 150, 300 GeV in the case of Ues = 0.2 with = lO^^GeV. 
In fig.^ the branching ratio is shown for Ues = 0.05 with = lO^^GeV. The Mr dependence is 
the same as the case of the quasi-degenerate neutrino masses. The predictions almost exceed the 
experimental bound as far as Ue 3 > 0.05, tan /3 > 10 and ~ lO^^GeV. This result is based on 
the assumption Mj^ = M^l, however, it is not guaranteed in the case of the inverse-hierarchical 
neutrino masses. 

For comparison, we show the branching ratio in the case of the hierarchical neutrino masses 
in fig.|^ It is similar to the case of the inverse-hierarchical neutrino masses. The branching 
ratio in the case of the degenerate neutrino masses is 10^ times smaller than the one in the 
inverse-hierarchical and the hierarchical neutrino spectra. 

The analyses in the previous section depend on the assumption of M^i = M112 = = Mr. 

In the case of the quasi-degenarate neutrino masses in eq. this complete degeneracy of 
may be deviated in the following magnitude without fine-tuning: 


Mk 


~ 1 ± 


m?. 


By using eq.(2.2), we obtain 


kk 

Ki 


~ 1 ± 


AmQ 

ml 


(3.6) 


(Yt Y,)2 i = (Yt Y, 




M- 


R-ocl 


+ A(YtY,) 


21 


(3.7) 


where the first term is the (YlYjy) 2 i element in eq.(^.4|), which corresponds to the while 

the second term stands for the deviation from it as follows: 

A/T ^ 

A(YJYj,)21 = —^ U2iUli^Jmui^Jmuj{e2R2iR2j + ■ (3.8) 

id 


In order to estimate the second term, we use £2 = 0.0001 and 63 = 0.01 taking account of 
62 — AmQ/2m^ and 63 ~ Arn^tm/^^^^) where mu = 0.3eV is put. Since mui — myj and 
Rij < 1 , we get 

A(Y);Y^)2i < 3.5 X 10“^ . (3.9) 


Taking this maximal value, we can estimate the branching ratio as follows: 

i?/2(non-degenerate Mfl) /2.6-|-3.5\^ 

B/;{degenerate Af.) 

Therefore, the enhancement due to the second term is at most factor 5. This conclusion does 
not depend on the specific form of R 

Consider the case of the inverse-hierarchical type of neutrino masses. We take 62 ~ 0.01 
with the similar argument of the quasi-degenerate type neutrino masses, because mui and mu 2 
are almost degenerate and 62 — AmQ/2Am^^j^ in this case. Then, we get 

A(y1;Y^)2i < 0.063 X 10-2 , 


(3.11) 












where we assume £2 > £3 and use 771^,3 ~ 0, rriyi ~ my 2 — 0.054eV and Rij < 1. Taking the 
maximal value, we get 

5/2(non-degenerate Mji) / 2.7 +0.063 ^ 
iifl(degenerate Jf.) “ I, 2.7 j " 

Thus, the effect of the A(YJYj,) 21 is very small in the case of the inverse hierarchical neutrino 
masses. 
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Degenerate type MR=10*‘^GeV 



Figure 1: BR(/i ^ e + 7 ) in the case of the 
degenerate neutrino masses. 



Figure 2: BR(// —> e + 7 ) in the case of the case 
of degenerate neutrino masses. 


Inverse-hierarchical type MR=10*^GeV 



Figure 3: BR(;U ^ e + 7 ) in the case of the 
inverse-hierarchical neutrino masses. 


Inverse-hierarchical type MR=10''‘GeV , Ue3=0.05 



Figure 4: BR(/i —> e + 7 ) in the case of the 
inverse-hierarchical neutrino masses. 


Hierarchical type MR=10'‘^GeV 



Figure 5: BR(/i —> e -|- 7 ) in the case of the 
hierarchical neutrino masses. 
























































